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ABSTRACT 


An  objective  layer  of  maximum  wind  (LRMW)  analysis  technique  is  described 
and  evaluated.  The  evaluation  indicates  that  the  technique  represents  a  significant 
improvement  over  the  previously  developed  level  of  maximum  wind  (LMW)  analysis 
technique  [Spiegler,  D.  B.,  and  J.  T.  Ball,  et  al.,  1965:  Techniques  for  objective 
hemispheric  analysis  and  prediction  of  the  jet  stream.  Technical  Report  7463-176 
(ESD-TR-65-13),  The  Travelers  Research  Center,  Inc.]. 

A  categorization  procedure  for  wind  profiles  that  results  in  nine  jet  stream 
categories  is  designed  and  used  for  the  derivation  of  regression  equations  that  supply 
the  initial -guess  fields  for  each  of  five  LRMW  parameters.  The  stratification  of  the 
profiles  into  nine  categories  is  a  refinement  of  the  previous  seven  categories  developed 
for  the  LMW  analysis  technique.  The  initial-guess  LRMW  equations  are  stable  in 
tests  with  independent  data  and  are  capable  of  specifying,  very  well,  the  general 
characteristics  of  the  LRMW  profile.  They  also  provide  realistic  values  for  the 
LRMW  parameters  over  no-data  areas  that  are  consistent  with  the  entire  analysis  area. 

The  analysis  technique  locates  jet  cores  between  grid  points  and  generates 
observations  along  these  cores  by  using  models  of  horizontal  jet  profiles. 

Thickness  of  the  LRMW  as  one  of  the  analysis  fields  in  conjunction  with  the  wind- 
speed  maximum,  mean  height  of  the  LRMW,  and  vector  shears  below  and  above  the 
LRMW  enable  the  generation  of  a  three-dimensional  picture  of  the  wind  field  in  the 
vicinity  of  the  jet  stream. 

Over  the  data  areas,  where  a  comparison  between  objective  and  subjective 
analyses  can  be  made,  the  objective  LRMW  analyses  compare  favorably  with  the 
subjective  analyses. 
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SECTION  I 


INTRODUCTION 

An  objective  hemispheric  analysis  technique  for  five  level-of-maximum-wind 
(LMW)  parameters  was  described  and  evaluated  in  a  previous  report  [6].  The  five  LMW 
parameters  are:  the  wind  direction  at  the  LMW,  the  wind  speed  at  the  LMW,  the  height 
of  the  LMW,  and  the  scalar  shears  10,000  ft  above  and  below  the  LMW.  The  evaluation 
indicated  that  the  analysis  technique  produces  high  quality  analyses,  particularly  of  the 
wind  speed,  but  that  there  are  a  few  areas  of  the  hemisphere  near  jet  streams  where  the 
height  of  the  LMW  and  the  vertical  shears  are  not  specified  properly.  These  areas  are 
usually  associated  with  wind  profiles  that  exhibit  a  thick  layer  of  high  wind  speeds  (i.e., 
cases  where  there  is  more  than  one  distinct  maximum,  such  as  in  areas  where  both  the 
polar  and  subtropical  jet  are  observed  Because  the  analysis  technique  requires  a  level 
of  maximum  wind,  a  particular  level  is  selected  somewhere  in  the  thick  layer,  although 
the  level  selected  is  not  meaningful  for  these  cases.  The  analyses  and  verifications  of 
the  height  and  vertical  shears  are  sensitive  to  where  in  the  layer  the  LMW  is  specified. 
Because  of  this,  it  was  concluded  that  the  analysis  technique  should  be  refined  to  pre¬ 
sent  a  more  realistic  three-dimensional  picture  of  the  jet  stream  and  its  associated 
parameters.  To  accomplish  this  refinement,  a  layer  of  maximum  wind  (LRMW)  was 
defined  (similar  to  that  discussed  by  Reiter  [5])  and  the  analysis  of  six  LRMW  parameters 
was  performed.  These  parameters  are: 

(a)  the  mean  altitude  of  the  LRMW  (bounded  by  levels  where 
the  wind  speed  is  85%  of  the  maximum  speed  in  the  layer),  Z  (LR), 

(b)  the  thickness  of  the  LRMW,  h  (LR), 

(c)  the  maximum  wind  speed  in  the  LRMW,  Wg^ 

(d)  the  mean  wind  direction  in  the  LRMW,  V^(LR), 

(e)  the  vertical  vector  shear  below  the  LRMW,  and 

(f)  the  vertical  vector  shear  above  the  LRMW,  S  . 

The  mean  altitude  of  the  layer  of  maximum  wind  replaces  the  height  of  the  level 
of  maximum  wind  in  the  analysis  technique.  The  thickness  of  the  LRMW  is  a  new  field. 
Vector  shears  above  and  below  the  LRMW7,  instead  of  scalar  shears  above  and  below  the 
LMW,  complete  the  more  realistic  three-dimensional  analysis.  The  LRMW  analysis 
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technique  is  similar  to  that  described  for  the  LMW  [6];  that  is  (a)  initial-guess  fields 
are  supplied  by  using  an  established  relationship  between  the  LRMW  parameters  and 
the  information  available  from  constant-pressure-surface  analyses  of  wind,  temperature, 
and  geopotential  height,  and  (b)  the  initial-guess  fields  are  adjusted  using  station  values 
of  the  LRMW  parameters  (computed  by  a  preprocessor  program)  and  the  successive 
approximation  analysis  technique  (SAT)  [1]. 

Section  II  describes  the  experiments  performed  and  the  analysis  technique,  and 
Section  III  discusses  the  results  of  the  tests  and  the  evaluation  of  the  technique. 
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SECTION  II 


EXPERIMENTS  FOR  LAYER  OF  MAXIMUM  WIND 
INITIAL-GUESS  EQUATIONS.  AND  THE  PROCEDURE  FOR  ANALYSIS 

1.  Selection  of  LRMW  Categories 

In  previous  work  on  the  objective  analysis  of  the  jet  stream  [6],  seven  jet-stream 

categories  were  defined  based  on  wind  speeds  at  200  and  300  mb  [W  (2)  and  W  (3)], 

s  s 

latitude,  and  a  combination  of  both  wind  speeds  and  latitude.  These  categories  were: 
no  jet,  subtropical  jet'*’,  subtropical— polar  jet'*’,  and  polar  jet’*'. 

The  “subtropical— polar”  jets  are  characterized  by  either  a  single  maximum 
somewhere  between  200  and  300  mb,  or  by  multiple  maxima  (e.g.,  polar  jet  in  the 
vicinity  of  300  mb  and  subtropical  jet  in  the  vicinity  of  200  mb).  Thus,  it  was  advanta¬ 
geous  to  further  separate  the  two  subtropical-polar  categories  into  single-  and  multiple- 
type  wind  maxima.  The  resulting  nine  categories  and  the  associated  criteria  for  each 
are  listed  in  Table  I. 

Latitude,  which  had  previously  been  included  as  one  of  the  criteria,  was  omitted 
because  it  was  noted  that  the  jets  to  the  north  of  strong  mid-  to  high-latitude  ridges, 
or  blocking  anticyclones,  had  characteristics  of  subtropical  jets,  even  though  located 
north  of  45°  (e.g.,  the  height  of  the  maximum  wind  speed  was  around  40,000  ft).  Con¬ 
versely,  the  polar  jet  was  sometimes  observed  in  the  subtropics  at  a  latitude  below 
30°  (30°N  and  45°N  had  previously  been  the  limiting  latitudes  for  polar  and  subtropical 
jets  [6]). 

2 

Thickness  of  the  LRMW  ,  h(LR),  was  included  in  the  criteria  for  the  subtropical- 
polar  jets  for  the  reasons  previously  discussed.  A  12,000-ft  thickness  was  chosen  as 
the  separation  limit  between  single  and  multiple  maxima  after  considerable  testing 
with  8,000-  and  10,000-ft  thicknesses. 

Criteria  established  for  the  nine  jet-stream  categories  should  help  to  overcome 
the  problem  of  accurately  specifying  the  shears  and  mean  height  of  the  layer  of  maxi¬ 
mum  wind  caused  by  a  thick  LRMW. 

*Two  categories:  “weak— moderate”,  and  “strong”. 

2 

Thickness  of  the  LRMW  is  defined  by  the  height  difference  between  the  two  levels, 
above  and  below  the  wind  speed  maximum,  at  which  the  wind  speed  is  equal  to  85%  of 
the  wind  speed  maximum. 
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2. _ LRMW  Screening  Regression  Experiments 

The  initial-guess  fields  to  the  LRMW  analyses  are  provided  by  regression  equations 
derived  for  the  nine  jet-stream  categories  described  in  Section  1. 

To  derive  these  regression  equations,  a  three-month  sample  of  winter  upper-air 
data  (Dec  1963,  Jan  and  Feb  1964)  received  from  the  National  Meteorological  Center 
(NMC)  was  completely  error  checked  and  processed  through  three  programs  before  the 
screening  experiments  were  run.  The  first  of  these  programs  selects  and  computes 
LRMW  parameters  from  rawinsonde  and  radiosonde  station  data.  The  second  processor 
program  computes  a  list  of  specifiers,  which  may  be  selected  for  use  in  the  derived 
LRMW  equations.  The  specifiers  are: 

(a)  temperature,  winds,  and  heights  at  mandatory  constant-pressure 
surfaces  from  700  to  100  mb; 

(b)  thickness,  temperature  differences,  shears,  and  mean  temperatures 
of  layers  bounded  by  mandatory  constant -pressure  surfaces; 

(c)  thicknesses  between  the  mean  height  of  the  LRMW  and  the  height 
of  constant -pres  sure  surfaces  from  500  to  100  mb; 

(d)  wind  shears  between  the  maximum  wind  speed  and  the  wind  speed 
at  mandatory  constant-pressure  surfaces;  and 

(e)  an  estimated  thickness  of  the  LRMW  (computed  from  constant- 
pressure-surface  winds  and  heights;  the  details  of  the  computation  are  in 
Section  3). 

These  specifiers  are  listed  in  Table  II. 
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TABLE  II 

SPECIFIERS  USED  IN  LRMW  SCREENING  REGRESSION  EXPERIMENTS* 


Z(7),  T(7),  Wd(7),  Wg(7) 

Z(5),  T(5),  Wd(5),  Wg(5) 

Z(4),  T(4),  Wd(4),  Wg(4) 

Z(3),  T(3),  Wd(3),  Wg(3) 

Z(l),  T(l),  Wd(l),  Wg(l) 

h(7— 5),  AT(7— 5),  AW  (7-5),  AT/An(7-5) 

s 

h(5— 4),  AT(5— 4),  W  (5-4),  AT/An(5-4) 

s 

h(5— 3),  AT(5— 3),  AW  (5-3),  AT/An(5-3) 

s 

h(4— 3),  AT(4— 3),  AW  (4-3),  AT/An(4-3) 

s 

h(2— 1),  AT(2— 1),  AW  (2-1),  AT/An(2-1) 

s 

h(5 — M),  W  (5— M),  S 
s  o 

h(3—  M),  W  (3— M),  S 

S  j 

h(l-M),  W  (1— M),  S, 
s  i 

W  (L) 
s 

h(LR),  he(LR),  Wd(LR),  Z(LR),  Sfc,  ?  .  S^, 

*Z  =  altitude;  2  -  mean  altitude;  T  =  temperature;  AT  =  temperature 
difference;  W^  =  wind  direction;  Wg  =  wind  speed;  AWS  =  wind  shear; 
h  =  thickness;  he  =  estimated  thickness;  AT/ An  =  thermal  wind; 

Sx  =  [Ws(x-M)]/[Z(x)-Z(LR)];  Sa  and  Sjj  =  vector  shears  above  and  below; 

and  =  maximum  Sa  and  S^;  parenthetic  numbers  refer  to  level  in 
100  mb. 
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The  third  processor  program  separates  the  data  into  the  nine  categories  discussed 
in  the  previous  section. 

3.  The  LRMW  Analysis  Procedure 

The  regression  equations  for  the  nine  jet-stream  categories  provide  the  initial 
guesses  (grid-point  values)  for  the  analyses  of  all  LRMW  fields  except  wind  direction. 
The  procedures  for  obtaining  grid-point  values  for  each  of  the  LRMW  parameters  are 
as  follows. 

3.1  Wind  direction  for  the  LRMW.  W^(LR) 

(a)  For  subtropical-jet  categories  (2  and 3), 

Wd(LR)  =  Wd(2) 

because  height  of  the  subtropical  jet  is  generally  near  200  mb. 

(b)  For  polar-jet  categories  (8  and  9), 

Wd(LR)  =  Wd(3) 

because  height  of  the  polar  jet  is  generally  near  300  mb. 

(c)  For  subtropical — polar  and  no-jet  categories  (4—7  and  1), 

Wd(LR)  =  Wd(3)  +  Wd(2) 

2 

because  a  subtropical— polar  jet  is  generally  located  between  300  and  200  mb; 
for  no- jet  categories,  “maxima”  may  be  anywhere  in  the  sounding,  and  the  average 
of  wind  directions  at  300  and  200  mb  is  a  good  approximation. 

3.2  Maximum  wind  speed  within  the  LRMW,  W 

The  jet-stream  category  for  each  grid  point  is  determined  in  order  to  generate 
initial  guesses  for  the  maximum  wind  field.  This  is  not  a  straightforward  procedure 
for  categories  4  through  7,  because  thickness  of  the  LRMW  (h(LR))  is  part  of  the  cate¬ 
gory  criteria  and  also  an  LRMW  parameter  to  be  analyzed,  but  is  not  available  at  this 
point.  Categories  1,  2,  3,  8, and  9,  on  the  other  hand,  are  uniquely  determined  by  300- 
and  200-mb  wind  speeds  and  their  differences. 

Category  4  is  different  from  Category  6  only  in  the  thickness  criteria  (similary 
for  Categories  5  and  7).  The  solution  to  the  problem  is  the  computation  of  an  estimated 
thickness  of  the  LRMW  [h  (LR)J  to  classify  uniquely  grid  points  for  Categories  4  and  6 
and  Categories  5  and  7.  The  only  information  available  (at  all  grid  points)  for  the 
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computation  is  constant-pressure-surface  data.  The  estimated  thickness  is  computed 
as  follows. 

The  level  of  maximum  wind  (LMW)  is  assumed  to  be  at  36,000  ft  for  Categories 
4,  5,6,  and  7  (as  this  was  the  mean  height  to  the  nearest  thousand  feet  of  the  LMW  for 
subtropical— polar  jet  categories  {  1840  cases}  [6]).  For  a  grid  point  that  may  be  in 
either  Category  4  or  6,  based  on  the  300-  and  200-mb  wind  criteria,  the  wind-speed 
maximum  is  computed  for  both  categories  and  the  average  value  is  used  as  an  aid  in 
obtaining  an  estimated  thickness.  Eight-five  percent  of  the  averaged  wind-speed  maxi¬ 
mum  defines  the  boundaries  of  the  layer  of  maximum  wind.  A  simple  interpolation  be¬ 
tween  the  heights  of  constant-pressure  surfaces  above  and  below  the  LMW  and  their 
associated  wind  speeds  determines  the  heights  at  which  the  wind  speed  is  85%  of  the 
maximum  speed.  The  estimated  thickness  of  the  LRMW  is  obtained  by  subtracting  the 
interpolated  height  of  the  lower  boundary  of  the  layer  from  the  interpolated  height  of 
the  upper  boundary  of  the  layer.  The  proper  category  for  the  grid  point  can  then  be 
determined  with  the  estimated  thickness.  A  similar  procedure  is  employed  for  grid 
points  that  initially  could  be  in  Category  5  or  7. 

The  suitability  of  this  method  for  properly  categorizing  a  grid  point  was  tested  on 
the  3-month  upper-air  data  sample.  The  number  of  cases  placed  in  the  proper  category 
using  the  estimated  thickness  is  given  in  Table  III;  the  actual  number  of  cases  in  a 
category  was  determined  by  using  winds-aloft  station  data  (dependent)  to  arrive  at  the 
true  thickness  [h(LR)]. 

The  results  in  Table  III  indicate  that  the  method  of  estimating  the  thickness  of  the 
LRMW  works  well. 

With  all  grid  points  now  uniquely  defined  by  category,  the  initial  guess  of  W 

o  1V1 

is  generated  by  the  appropriate  regression  equations.  The  initial  guess  is  “corrected” 
by  using  the  processed  station  observations  and  the  successive  approximation,  analysis 
(SAT)  technique  [1]. 

After  a  specified  number  of  passes  through  the  data,  the  jet  cores  are  located 
(first  at  grid  points,  and  then  between  grid  points  through  a  parabola  fitting  routine). 
“Observations”  are  generated  along  core  segments  by  using  Endlich  and  McLean’s 
jet-stream  model  [2],  which  relates  the  wind  speed  at  the  core  to  the  wind  speed  at  some 
specified  perpendicular  distance  from  the  core.  This  model  was  first  derived  from 
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aircraft  flights  in  the  southeastern  U.S.  and  results  were  recently  corroborated  by 
similar  flights  into  jet  streams  in  the  central  U.S.  [3]. 

TABLE  HI 

NUMBER  OF  CASES  PLACED  IN  PROPER  CATEGORY 
USING  h  (LR)— DEPENDENT  DATA 


Jet-stream 

category* 

No.  of  cases 
hg(LR) 

computation 

No.  of  cases, 
h(LR) 

computation 

% 

correct 

4 

798 

935 

85.3 

5 

231 

267 

86.5 

6 

534 

779 

75.5 

7 

368 

468 

78.6 

Total 

1931 

2449 

78.8 

*See  Table  I  for  definitions. 

An  additional  pass  is  then  made  using  both  the  generated  data  and  the  actual 
station  data,  and  a  relatively  small  influence  radius.  The  generated  wind-speed 
“observations”  may  be  assigned  a  weight  relative  to  the  actual  data  for  this  pass. 

Verification  of  W 
- sM 

The  areal -mean  rms  (ARMS)  error  method  is  used  for  verification  purposes  [7]. 
Although  it  was  described  in  the  previous  report  [6 J  (and  elsewhere),  it  is  described  here 
because  verification  is  considered  one  of  the  important  aspects  of  the  analysis  evalua¬ 
tion.  The  ARMS  error  method  approximates  an  area  integral  of  the  error  over  the 
analysis  area  by  withholding  some  percentage  of  the  observations  and  verifying  them 
against  the  withheld-  and  analysis-station  observations,  which  are  equally  weighted.  The 
withheld-station  observations  approximate  the  maximum  error  points  in  the  analysis 
field,  but  the  number  of  withheld  stations  required  by  the  ARMS  error  method  usually 
does  not  equal  the  number  of  maximum  error  points.  Therefore,  the  withheld  stations 
were  selected  to  represent  equally  the  entire  spectrum  of  maximum  errors  over  the 
analysis  area.  The  assumption  is  that  the  maximum  error  is  inversely  proportional  to 
the  density  of  the  data. 
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The  “density”  of  each  station  is  computed  prior  to  the  analysis  (see  [7]  for  details 
of  the  density  computation),  and  withheld  stations  are  selected  on  the  basis  of  three 
maximum-wind-speed  and  three  station-density  categories  as  outlined  in  Table  IV. 


TABLE  IV 

SELECTION  OF  WITHHELD  STATIONS  ON  THE  BASIS  OF  STATION 
DENSITY  AND  WIND  SPEED  MAXIMUM 


No.  of  withheld  stations 

Density 

(P) 

A 

[W  >  100  knots] 

sM 

B 

[W  =  60-100  knots] 

sM 

C 

[W  <  60  knots] 

sM 

Low 

(0.11—0.32) 

3 

4 

3 

Medium 

(0.33-0.43) 

4 

3 

3 

High 

(£0.44) 

3 

3 

4 

3.3  Thickness  of  the  LRMW,  h(LR) 

This  field  is  analyzed  at  all  grid  points  except  those  in  Category  1.  The  initial 
guess  for  h(LR)  is  generated  by  regression  equations  after  determining  the  jet-stream 
category  for  each  grid  point.  The  range  of  thickness  for  Category-1  data  is  too  large 
to  produce  meaningful  values.  The  SAT  technique  is  then  used  with  the  station  data  for 
h(LR)  to  adjust  the  initial-guess  field.  No  station  is  selected  for  corrections  to  grid 
points  where  Wgj^  <  60  knots  at  the  station. 

Verification  of  h(LR) 

Thickness  of  the  LRMW  is  verified  using  the  ARMS  error  method  and  the  identical 

set  of  withheld  stations  used  in  the  verification  of  W  „  with  the  exception  of  verification 

sM 

Category  C  [W^  <60  knots]. 

3.4  Mean  height  of  the  LRMW,  Z(LR) 

The  procedure  for  generating  the  initial-guess  field  to  the  analysis,  the  “correction” 
to  grid  points,  and  the  verification  is  identical  to  that  for  the  thickness  of  the  LRMW. 
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3.5  Vector  shears  below  and  above  the  LRMW,  and 

The  vector  shears  below  and  above  the  LRMW  are  analyzed  separately,  with  the 
analysis  for  the  shear  below  performed  first. 

The  initial  guesses,  the  analysis,  and  the  verifications  are  the  same  as  those  for 
the  maximum  wind  speed  in  the  LRMW  (i.e.,  analysis  and  verifications  at  all  grid 
points). 

On  the  shear-below  and  shear-above  printout  charts,  station  locations  where 
maximum  vector  shears  (anywhere  in  the  wind  profile)  exceed  10  knots  per  thousand 
feet  are  indicated  by  code  letters.  The  code  letters  correspond  to  a  list  of  stations 
that  is  printed  with  the  corresponding  maximum  shear  values  and  the  levels  at  which 
the  maximum  shear  is  observed. 


SECTION  III 


RESULTS  FROM  LRMW  EXPERIMENTS 

The  quality  of  the  LRMW  analysis  technique  was  determined  from:  (a)  testing  the 
equations  that  provide  the  initial  guesses  to  the  analyses  with  independent  station  data, 
and  (b)  evaluating  a  series  of  objective  analyses  using  both  the  ARMS  error  method 
and  comparisons  with  subjective  analyses  of  LRMW  parameters.  The  test  and  evalua¬ 
tion  of  these  items  is  described  in  the  following  two  sections. 

4.  LRMW  Initial-guess  Equations— Test  and  Evaluation 

The  list  of  possible  specifiers  for  the  LRMW  parameters  includes  those  specifiers 
directly  available  on  a  routine  basis  from  data  tapes  (e.g.,  winds  and  temperatures  at 
constant -pres sure  surfaces),  those  that  are  computed  from  constant-pressure-surface 

wind  and  height  data  (e.g.,  thermal  winds,  thicknesses),  and  those  that  contain  one  or  two 

3 

LRMW  parameters  ,  either  alone  or  as  some  arithmetic  expression  with  constant- 
pressure-surface  data  (e.g.,  wind  speed  maximum,  or  shear  between  the  wind  at  the 
LMW  and  at  a  constant-pressure  surface). 

Several  screening  regression  experiments  were  performed  for  each  specificand 
in  each  category.  The  primary  objective  in  designing  several  experiments  was  to 
determine  what  combination  of  possible  specifiers  would  yield  the  best  results  (some 
experiments  used  only  directly  available  parameters  as  possible  specifiers;  others 
used  only  computed  parameters;  and  still  others  used  various  combinations  of  directly- 
available  and  computed  parameters  as  possible  specifiers). 

Table  V  lists  the  specificands,  selected  specifiers  and  their  associated  per  cent 
reduction  of  variance  (%  red.),  number  of  cases,  rms  errors,  and  standard  deviations 
for  each  jet  category  for  both  dependent  and  independent  data.  Figure  1  illustrates  the 
dependent  and  independent  rms  errors  for  each  LRMW  pararmeter.  A  study  of  the 
table  reveals  that: 


3 

These  specifiers  assume  (for  a  particular  LRMW’  parameter)  that  the  analysis  of 
the  parameters  proceeds  in  a  particular  sequential  order  (which  it  does— sec  Section  3). 
Thus,  the  wind  speed  may  be  a  possible  specifier  for  all  other  LRMW'  parameters  be¬ 
cause  it  is  analyzed  first. 
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TABLE  V  (Continued) 
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 s  s 

Std.  dev. 

Ind. 

14.7 

2784 

6094 

2.39 

2.78 

18.0 

3459 

4983 

3.64 

3.28 

12.6 

3520 

5341 

2.17 

2.19 

Dep. 

15.9 

2744 

5516 

2.60 

2.75 

20.3 

2706 

4211 

3.15 

3.77 

12.5 

3311 

4583 

2.22 

2.13 

Rms  error 

1 ! 

Ind. 

9.5 

2383 

5089 

2.05 

2.42 

11.3 

2730 
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1 
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TABLE  V  (Continued) 
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Fig.  1.  Root-mean-square  errors  of  LRMW  parameters. 
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Subtropical -polar 
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Fig.  1.  Continued. 
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(a)  The  equations  are  stable  with  independent  data  because 
independent  rms  errors  are  frequently  near  the  rms  dependent  errors; 
where  the  independent  errors  are  higher  (or  lower)  than  the  dependent 
errors,  the  standard  deviation  of  the  independent  sample  is  also  higher 

(or  lower),  except  for  S  ,  Category  6. 

a  4 

(b)  In  general,  the  specifiers  selected  by  the  screening  program 

are  those  that  one  might  expect  to  be  selected  from  synoptic  experience 
with  vertical  wind  profiles,  i.e.,  (1)  the  specifiers  selected  for  the  maxi¬ 
mum  wind  speed  for  all  categories  are  one  or  more  of  the  directly  avail¬ 
able  constant -pres sure -surface  winds  at  300,  200,  and  100  mb.  Consider¬ 
able  skill  is  demonstrated  in  specifying  the  maximum  wind,  as  evidenced 
by  the  large  %  red.in  both  the  independent  and  dependent  samples; 

(2)  for  thickness  of  the  LRMW,  the  estimated  thickness  is  one  of  the 
selected  specifiers  for  seven  of  eight  LRMW  categories  (thickness  not 
specified  for  no  jet,  Category  1);  (3)  for  vector  wind  shears  below  and 
above  the  LRMW,  the  computed  wind  shears  between  constant-pressure 
surfaces  and  between  the  maximum  wind  level  and  constant-pressure 
surfaces  are  invariably  among  the  selected  specifiers;  (4)  the  mean  height 
of  the  LRMW  is  best  specified  by  a  combination  of  wind  and  temperature 
information  at  500,  200,  and  100  mb  for  most  of  the  LRMW?  categories. 

(c)  There  is  a  significant  difference  in  the  mean  thicknesses  of 
the  LRMW  for  the  single-  and  multiple -maxima  subtropical— polar  jet 
categories,  i.e.,  the  mean  thicknesses  are  not  just  below  and  above  the 
limiting  value  of  12,000  ft  given  by  the  definitions;  the  mean  thicknesses  for 
the  multiple  jet  categories  (6  and  7)  are  above  twice  that  for  the  single  jet 

4 

The  cut-off  criteria  for  selection  is  a  combination  of  the  modified  F-test  described 
by  Miller  [4]  and  practical  considerations  involving  the  optimum  number  of  specifiers 
to  use  on  an  operational  basis  (e.g.,  a  specifier  may  have  been  computed  to  be  sig¬ 
nificant  by  the  modified  F-test,  but  may  not  have  been  included  with  those  listed  in 
Table  V  if  the  %  red.  and  the  lowering  of  the  dependent  rms  error  was  considered  to 
be  small  enough  to  make  little  difference  in  specification  accuracy,  considering  the 
increased  storage  and  computer  time  required  for  its  use). 
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categories  (4  and  5).  Thus,  the  stratification  into  single  and  multiple  jet 
categories  is  justified  by  the  data  and  is  a  desirable  revision  of  the  sub¬ 
tropical-polar  jet  categories  defined  in  the  earlier  report  [6],  where 
thickness  was  not  a  criterion. 

Table  VI  shows  the  percentage  of  cases  within  specified  error  limits  for  the 
independent-data  test  of  the  LRMW  equations.  The  error  limits  in  the  table  were 
selected  as  representing  a  range  of  errors  that  is  quite  low  considering  that  the  LRMW 
equations  were  derived  to  provide  initial-guess  fields  to  analyses  of  LRMW  parameters. 
It  can  be  seen  from  the  table  that  a  very  high  percentage  of  the  specification  errors 
fall  within  these  limits. 

4 .1 _ Subtropical  and  Polar  Jets  with  h  (LR)  >  12,000  ft 

It  was  recognized  that  the  thickness  of  the  layer  of  maximum  wind  could  exceed 
12,000  ft  for  a  wind  profile  that  clearly  indicated  either  a  subtropical  jet  or  a  polar  jet. 


TABLE  VI 

PERCENTAGE  OF  CASES  WITHIN  SPECIFIED  ERROR  LIMITS 
FOR  LRMW  PARAMETERS  (Independent  data) 


Jet-stream 

category 

Total 

number  of 

cases 

Error  Limits 

w 

sM 

±  10  knots 

h(LR) 

±  4000  ft 

Z(LR) 

±  4000  ft 

S 

b 

±  2  knots 

X‘ 

i 

533 

81.3 

— 

— 

77.8 

76.0 

2 

206 

88.8 

79.1 

89.3 

79.6 

72.8 

3 

97 

94.9 

87.6 

99.0 

80.3 

65.9 

4 

196 

74.0 

90.3 

73.5 

74.0 

66.9 

5 

49 

61.2 

87.8 

81.6 

61.2 

51.0 

6 

192 

93.7 

78.6 

79.1 

75.5 

62.1 

7 

95 

96.8 

85.3 

92.6 

70.5 

77.9 

8 

108 

81.5 

79.6 

91.6 

63.9 

63.0 

9 

28 

85.7 

96.4 

100.0 

53.5 

67.9 
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To  determine  the  frequency  of  occurrence  of  these  events,  the  dependent  data 
sample  for  Categories  2  and  3  (subtropical)  and  Categories  8  and  9  (polar)  were 
separated  into  cases  where  h(LR)  s  12,000  ft  and  cases  where  h(LR)  >  12,000  ft. 
Table  VII  shows  the  breakdown  of  cases. 

TABLE  VII 

NUMBER  AND  PERCENTAGE  OF  CASES  WHERE  h(LR)  >  12,000  ft 
FOR  JET-STREAM  CATEGORIES  2,  3,  8,  AND  9 
(Dependent  data) 


Jet-stream 

category* 

No.  of  cases 

Total  no. 
of  cases 

%  cases 

h(LR)  >  12,000  ft 

h(LR)  £  12,000  ft 

h(LR)  >  12,000  ft 

2 

741 

213 

954 

22.3 

3 

541 

217 

758 

28.6 

8 

687 

126 

813 

15.5 

9 

272 

43 

315 

13.7 

*See  Table  I  for  definitions. 


It  is  readily  apparent  that  thicknesses  greater  than  12,000  ft  occur  relatively  infre¬ 
quently  for  these  categories.  To  determine  if  the  equations  derived  for  the  entire 
category  would  perform  satisfactorily  on  the  subset  of  data  given  by  the  cases  where 
h(LR)  >  12,000  ft,  the  equations  were  tested  with  the  subset.  The  statistics  appear  in 
Table  VIII. 

An  examination  of  Table  VIII  and  a  comparison  of  it  with  Table  V  indicates  that 
the  rms  errors  are  very  close  to  or  lower  than  the  errors  for  all  cases  in  the  category 
for  all  LRMW  parameters  except  the  thickness  of  the  LRMW.  It  is  not  surprising  that 
the  thickness  is  not  specified  well  for  these  cases,  because  the  mean  thicknesses  are 
about  5000  to  nearly  7000  ft  higher  than  the  mean  for  the  total  number  of  cases  in  the 
categories. 

Considering  the  stability  of  the  equations  on  the  subset  of  data  for  four  of  the 
five  LRMW  parameters,  the  relatively  infrequent  occurrence  of  thickness  >  12,000  ft 
for  the  polar  and  subtropical  jet  categories,  and  the  resultant  small  data  sample  avail¬ 
able  from  one  entire  winter’s  data,  it  was  concluded  that  separate  equations  were  not 
warranted  for  these  cases. 
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TABLE  VIII 

RESULTS  OF  TESTING  SUBTROPICAL-  AND  POLAR- 
JET-STREAM  EQUATIONS  ON  CASES  WHERE  h(LR)  >  12,000  ft 

(Dependent  data) 


Jet-stream 

category 

Specificand 

No.  of 
cases* 

Mean 

Rms 

error 

Units 

w  ™ 
sM 

84.7 

5.8 

knots 

h(LR) 

15133 

4349 

ft 

2 

Z(LR) 

212 

41727 

2174 

ft 

*b 

3.96 

1.62 

knots  10-3  ft--*- 

S 

a 

-3.82 

1.92 

knots  10  3  ft  1 

W  TV/I 

sM 

135.8 

5.7 

knots 

h(LR) 

14723 

3390 

ft 

3 

Z(LR) 

216 

40308 

1554 

ft 

*b 

5.24 

1.77 

knots  10“3  ft--*- 

S 

a 

-5.88 

2.03 

knots  10  3  ft  1 

W  » ir 

sM 

84.2 

5.7 

knots 

h(LR) 

14795 

5150 

ft 

8 

Z(LR) 

126 

28744 

2885 

ft 

?b 

4.73 

2.43 

knot  10-3  ft-1 

s' 

a 

-4.25 

2.42 

,  .  -3  -1 

knots  10  ft 

W 

sM 

121.5 

6.6 

knots 

h(LR) 

14785 

5247 

ft 

9 

Z(LR) 

43 

30976 

2758 

ft 

?b 

6.66 

2.28 

knots  10-3  ft-l 

"s 

a 

-5.45 

2.08 

-3  -1 

knots  10  ft 

*Where  the  number  of  cases  differs  fron  that  given  in  Table  VII,  a  case  was 
eliminated  due  to  data  error. 
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4.2  Estimated  h(LR)  Resulting  in  Placing  Cases  in  Wrong  Category 

For  jet-stream  Categories  4  through  7,  part  of  the  category  definition  is  given 
by  the  thickness  of  the  LRMW.  For  reasons  discussed  in  Section  3,  it  is  necessary 
to  compute  an  estimated  thickness  of  the  LRMW.  Table  III  lists  the  number  of  cases 
placed  in  the  proper  category  using  h^LR),  where  the  total  number  of  cases  in  a 
category  was  determined  by  using  winds-aloft  station  data  (dependent)  to  arrive  at 
the  real  thickness.  About  79%  of  the  cases  were  in  the  proper  category  for  the  four  jet- 
stream  categories.  A  similar  determination  was  made  for  the  station  data  from  the 
independent  sample  (see  Table  IX).  Then,  the  equation  for  the  categories  into  which 

TABLE  IX 

NUMBER  OF  CASES  PLACED  IN  PROPER  CATEGORY  USING 
h  (LR)-INDE PENDENT  DATA 


Jet-stream 

category 

No.  of  cases, 
he(LR) 

computation 

No.  of  cases, 
h(LR) 

computation 

% 

correct 

4 

163* 

196 

83.2 

5 

39| 

49 

79.6 

6 

160J 

193 

82.9 

7 

8511 

95 

89.5 

Total 

447 

533 

83.9 

*33  cases  fell  into  Category  6. 
tlO  cases  fell  into  Category  7. 

1 33  cases  fell  into  Category  4. 

H10  cases  fell  into  Category  5. 

cases  were  misplaced  [using  he(LR)]  were  used  with  the  data  from  the  categories  into 
which  the  cases  should  have  been  placed  [according  to  the  use  of  h(LR)]  to  determine 
how  well  the  equations  would  perform  for  these  misplaced  cases.  The  results  are  given 
in  Table  X. 

It  can  be  seen  that  applying  the  LRMW  equations  to  the  cases  where  erroneous 
categories  are  assigned  leads  to  large  errors  in  the  estimation  of  Z(LR).  However,  the 
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TABLE  X 

ERRORS  RESULTING  FROM  ASSIGNING 
WRONG  CATE GORIE S— INDE PE NDE NT  DATA 


Jet-stream  category 

Specific  and 

No.  of 

Mean 

Rms 

Units 

Equation 

Data 

cases 

error 

w  M 
sM 

83.3 

3.3 

knots 

h(LR) 

9012 

9583 

ft 

6 

4 

Z(LR) 

33 

35848 

6208 

ft 

?b 

3.43 

2.42 

knots  10-^  ft--'- 

S 

a 

-3.06 

3.14 

knots  10  **  ft 

W 

sM 

114.9 

3.0 

knots 

h(LR) 

8669 

8724 

ft 

7 

5 

Z(LR) 

10 

35571 

4370 

ft 

?b 

2.99 

1.91 

knots  10-3  ft-l 

S 

a 

-4.20 

1.87 

,  "3  ,  -1 

knots  10  ft 

W  V, 

sM 

85.7 

6.1 

knots 

h(LR) 

15309 

8097 

ft 

4 

6 

Z(LR) 

33 

37871 

4211 

ft 

3.78 

1.99 

knots  10-3  ft-1 

S 

a 

-3.92 

2.38 

knots  10  ^  ft  ^ 

w  M 
sM 

126.1 

8.2 

knots 

h(LR) 

14809 

6849 

ft 

5 

7 

Z(LR) 

10 

37080 

4326 

ft 

*b 

5.12 

2.35 

knots  10~3  ft-1 

s' 

a 

-5.64 

1.38 

-3  -1 

knots  10  ft 
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other  LRMW  parameter  specifications  are  reasonably  good.  These  results  are  not  diffi¬ 
cult  to  interpret.  As  given  by  the  definitions  of  the  categories  and  the  means  of  the  LRMW 
parameters,  the  major  difference  between  Category  4  and  Category  6  is  the  thickness  of 
the  LRMW,  and  similarly  for  Category  5  and  Category  7.  Other  differences  in  LRMW 
parameters  for  these  categories  are  not  readily  detectable.  Thus,  when  the  mean  thick¬ 
ness  of  the  LRMW  is  about  8,000  ft  and  a  thickness-specification  equation  from  a  cate¬ 
gory  where  the  mean  thickness  is  about  16,000  ft  is  used,  a  large  specification  error  is 
likely. 

As  noted  in  Section  3,  computation  of  an  estimated  thickness  of  the  LRMW  must 
use  constant-pressure-surface  information.  The  large  number  of  cases  in  each  of 
Categories  4  through  7  warrants  the  use  of  separate  LRMW  equations  for  each  of  these 
categories.  The  small  percentage  of  times  ( <  2 0% )  that  erroneous  categories  are  assigned 
should  not  seriously  affect  the  initial  guess  for  four  of  the  five  LRMW  analyses,  but 
will  affect  the  initial-guess  field  for  thickness  of  the  LRMW  over  10  to  25%  of  the  grid 
points  where  the  300-  or  200-mb  wind  speed  is  ^  60  knots.  However,  the  analysis  will 
reduce  considerably  the  number  of  grid  points  affected,  because  the  station  data  will 
adjust  the  grid-point  values . 

4.3  Category  Modeling  of  Wind  Profiles  (Categories  4—7) 

It  is  well  known  that  an  almost  infinite  variety  of  vertical  wind  profiles  is  pos¬ 
sible.  A  categorization  of  wind  profiles  for  the  purpose  of  an  objective  hemispheric 
analysis  technique  for  LRMW  parameters  is,  of  necessity,  based  on  constant-pressure  - 
surface  data  which  limits  the  sophistication  of  the  categorization  procedure.  Our 
attempt  to  achieve  modeling  sophistication  was  made  by  introducing  thickness  of  the 
LRMW  as  one  of  the  criteria  for  defining  the  profiles  that  were  other  than  “pure” 
subtropical  or  polar  jets. 

The  logic  for  using  300-  to  200-mb  wind-speed  differences  as  another  part  of 
the  category  criteria  is  straightforward.  If  there  is  a  relatively  small  difference, 
either  the  maximum  wind  is  between  300  and  200  mb  or  there  are  maxima  near  both 
of  these  levels.  The  thickness  of  the  LRMW  criteria  is  designed  to  separate  these 
either/or  cases;  i.e.,  if  the  thickness  is  ^  12,000  ft,  the  maxima  should  be  between 
300  and  200  mb,  and  if  it  is  >  12,000  ft,  there  may  be  maxima  near  300  and  200  mb 


27 


with  a  minimum  in  between  or,  if  the  maximum  is  between  300  and  200  mb,  it  must  be 
a  broad  maximum  (blunt  nose  jet)  by  definition. 

Where  there  is  a  relatively  large  difference  between  the  300  and  200  mb  wind 
speeds  (Categories  2,  3,  8  and  9),  the  maximum  wind  is  probably  near  the  eonstant- 
pressure  surfaee  having  the  higher  of  the  wind  speeds. 

To  determine  how  closely  the  wind  profiles  from  Categories  4  through  7 
followed  the  logic  of  the  criteria  established  for  them,  twenty-five  profiles  were  plotted 
from  station  data  of  winds  aloft  and  eonstant-pressure-surfaee  winds. 

Table  XI  lists  the  number  of  profiles  that  fell  within  the  criteria  limits 
described  above  for  each  of  Categories  4  through  7. 


TABLE  XI 

LEVEL  OF  MAXIMUM  WIND  LOCATIONS  FOR 
CATEGORIES  4  THROUGH  7 


Jet- 

stream 

category 

No.  of  eases  Wg^j  located: 

Minimum  between 
300  and  200  mb 
(between  2  maxima) 

Below 
300  mb 

At 

300  mb 

Between 

300  and  200  mb 

At 

200  mb 

Above 
200  mb 

4 

2 

— 

16 

2 

5 

— 

5 

— 

2 

17 

— 

6 

— 

6 

3 

— 

12 

2 

8 

5 

7 

1 

3 

12 

5 

4 

3 

For  the  single  jets,  65%  of  the  (25)  profiles  eontain  maximum  winds  between  300 
and  200  mb.  There  are  eases  where  the  difference  between  the  300  and  200  mb  wind 
speeds  is  not  large  enough  to  place  the  profiles  in  the  subtropical-  or  polar-jet  cate¬ 
gories,  but  where  the  speed  increases  above  200  mb  or  below  300  mb  to  result  in  the 
maximum  being  located  outside  the  300-  to  200-mb  layer. 

For  the  multiple  jets,  just  under  50%  of  the  profiles  eontain  maxima  between 
300  and  200  mb,  and  about  40%  contain  maxima  at  or  above  200  mb.  Sixteen  percent  of 
the  cases  exhibit  a  minimum  between  200  and  300  mb.  These  cases  arc  in  the  statistics 
twice,  with  maxima  at  or  near  200  or  300  mb.  Figures  2  through  5  show  examples 
of  profiles  for  eaeh  category  and  maximum  wind-speed  classification. 
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Station  72  253 
Station  72  662 
Station  72  311 


01  ‘I9A9X  B9S  9AOCJB  ;q8l9H 
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Fig.  3.  Examples  of  wind  profiles:  Category  5,  Single,  strong  subtropical— polar  jet. 
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Fig.  4.  Examples  of  wind  profiles:  Category  6,  multiple,  moderate  subtropical— polar  jet. 
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Fig.  5.  Examples  of  wind  profiles:  Category  7,  multiple,  strong  subtropical— polar  jet. 


In  general,  the  wind  profiles  follow  the  logic  for  the  category  criteria,  but,  as 
noted ,  there  are  instances  where  a  profile  for  a  particular  LRMW  category  deviates 
from  the  “model”  conceived  for  that  category. 

The  major  point  being  made  in  this  discussion  of  category  modeling  of  wind  pro¬ 
files  is  that  the  wind  profiles  cannot  be  adequately  stratified  to  include  all  types  of 
profiles,  at  least  when  using  only  the  constant-pressure-surface  information  (as 
required  by  the  nature  of  the  objective  LRMW  analysis  technique).  This  is  not  con¬ 
sidered  serious,  because  the  category  equations  supply  only  the  initial-guess  field 
for  the  LRMW  parameters,  and  the  initial  guess  is  “corrected”  over  the  data  areas  by 
the  analysis  procedure. 

It  is  concluded  that  the  LRMW  categories  defined  here  represent  the  near 
optimum  stratification  of  wind  profiles,  considering  that  the  definitions  must  be  based 
on  constant-pressure-surface  information. 

4.4  Verification  of  Different  Types  of  Wind  Profiles  (Categories  4—7) 

The  wind  profiles  in  Figs.  2  through  5  were  chosen  to  illustrate  some  of  the 
variety  of  wind  profiles  that  are  possible  within  particular  categories.  Subsequent  to 
their  selection  it  was  decided  to  look  for  the  individual  verifications  to  determine 
whether  the  equations  were  capable  of  “capturing”  some  of  the  differences  in  the  wind 
profiles.  For  each  of  the  profiles  in  Figs.  2  through  5,  it  is  possible  to  draw  portions 
of  the  wind  profiles  from  the  information  specified  by  the  LRMW  equations.  Figures 
6  through  9  represent  the  actual  observed  profiles  and  portions  of  the  profiles  con¬ 
structed  from  the  values  given  by  the  equations  for  the  LRMW  parameters.  On  these 
figures  the  maximum  wind  speed  in  the  LRMW  is  plotted  at  the  level  of  the  mean  height 
of  the  LRMW,  but  it  does  not  necessarily  follow  that  the  maximum  speed  is  actually  at 
that  level— it  may  be  anywhere  in  the  layer.  A  study  of  the  figures  and  specification 
errors  reveals  that  the  specification  of  LRMW  parameters  as  given  by  the  equations 
is,  in  general,  very  good.  Discussing  the  results  by  individual  categories: 

Category  4:  single, moderate  subtropical— polar  jet 

The  equations  are  seen  to  be  capable  of  specifying  the  mean  height  of  the  LRMW 
[Z(LR)]  above  200  mb,  below  300  mb,  and  inbetween  300  and  200  mb.  For  these  single 
jets,  it  may  be  assumed,  in  most  cases,  that  the  maximum  is  located  very  near  Z(LR) 
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Fig.  6.  LRMW  profiles  constructed  from  equations  and  observed  data— Category  4 
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Wind  speed,  knots 
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Wind  speed,  knots 
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Fig.  7.  LRMW  profiles  constructed  from  equations  and  observed  data — Category  5. 
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Wind  speed,  knots 


(c)  Station  72  311 
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Fig.  7.  Continued. 
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Fig.  8.  LRMW  profiles  constructed  from  equations  and  observed  data— Category  6  (boundaries  of  LRMW  connected 
by  solid  vertical  line  to  indicate  thickness). 
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Fig.  8.  Continued. 
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Fig.  8.  Continued. 
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Wind  speed,  knots 


o 

<D 


O 

CD 


O 

m 


o  o 

Tf  CO 


e 


01  ‘PA8X  B98  8AoqB 


o 

<N 


o 


44 


Fig.  9.  LRMW  profiles  constructed  from  equations  and  observed  data— Category  7  (boundaries  of  LRMW 
connected  by  solid  vertical  line  to  indicate  thickness). 
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Fig.  9.  Continued. 
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Wind  speed,  knots 


except  where  other  information  makes  this  unlikely  [e.g.,  Station  78  089  has  a  wind 
speed  of  86  knots  at  the  200-mb  level,  which  is  within  1000  ft  of  Z(LR);  thus,  it  is 
nearly  impossible  for  the  maximum  of  107  knots  to  be  at  the  specified  Z(LR)  of 
40557  ft.  With  the  information  of  the  boundaries  of  the  LRMW  given  by  the  thickness  in 
conjunction  with  the  known  200-mb  wind  speed,  it  is  very  likely  that  the  maximum  wind 
speed  is  located  between  200  mb  and  the  upper  boundary.  This  approximates  the  actual 
profile  very  well].  For  Station  47  582,  Z(LR)  is  specified  to  be  at  a  level  lower  than 
300  mb,  but  the  available  500-mb  wind  speed  (which  is  higher  than  the  maximum  speed 
specified  by  the  equations,  and  16  knots  greater  than  the  300-mb  speed)  strongly  sug¬ 
gests  that  the  maximum  is  closer  to  500  mb  than  300  mb, and  that  Z(LR)  is  probably 
specified  somewhat  high.  There  is  no  information  to  indicate  that  the  shears  are  in 
error;  thus,  the  profile  constructed  from  the  LRMW  equations  can  be  subjectively 
displaced  to  a  lower  elevation  somewhere  between  500  and  300  mb  to  give  a  better 
approximation  of  the  true  profile.  Station  72  655  exhibits  a  close  relationship  between 
the  observed  profile  and  that  specified  by  the  LRMW  equations. 

Category  5:  single,  strong  subtropical— polar  jet 

The  profiles  obtained  from  the  LRMW  information  approximate  rather  closely 
the  observed  profiles,  despite  the  observed  range  of  the  mean  height  of  the  LRMW  from 
29551  to  41213  ft.  Although  the  maximum  wind  speed  for  Station  72  311  was  not  speci¬ 
fied  to  be  as  high  as  observed,  most  of  the  sharpness  of  the  profile  was  indicated  by 
the  strong  vertical  shears  and  the  relatively  thin  layer  of  maximum  wind  given  by  the 
LRMW  equations. 

Category  6:  multiple,  moderate  subtropical— polar  jet 

In  this  category  the  vector  shears  above  and  below  the  LRMW,  as  well  as  the 
other  LRMW  parameters,  are  specified  quite  well  (see  Fig.  8).  The  maximum  will 
frequently  not  differ  by  an  appreciable  amount  from  the  300-  and  200-mb  wind  speeds; 
i.e.,  the  definition  of  the  category  precludes  a  sharpness  to  the  profile.  The  details 
within  the  thick  LRMW,  of  course,  are  not  specified. 


Category  7:  multiple,  strong  subtropical— polar  let 

The  comments  made  for  Category  6  apply  here.  LRMW  parameters  are  specified 
very  well. 

The  major  outcome  of  the  verification  results  for  these  assorted  profiles  is  an 
increase  in  confidence  in  the  ability  of  the  LRMW  equations  to  closely  approximate  the 
observed  values  of  the  LRMW  parameters  despite  the  incongruity  of  profile  types 
within  a  category. 

5.  LRMW  Objective  Analysis  Technique— Test  and  Evaluation 

The  analysis  procedure  described  in  Section  II  was  used  with  data  from  12  obser¬ 
vation  times  in  January  1963  (a  time  period  not  within  the  three  month  sample  used  to 
derive  the  initial-guess  equations).  The  LRMW  analyses  for  the  12  observation  times 
were  verified  using  the  ARMS  error  method.  The  evaluation  consisted  of  (a)  an  inter¬ 
pretation  of  the  verification  statistics  and  comparison  with  verification  statistics  from 
LMW  analyses  described  in  the  previous  report,  (b)  a  comparison  of  the  computer 
analyses  with  hand-analyzed  charts,  and  (c)  a  determination  of  the  consistency  between 
LRMW  analyses  for  each  observation  time. 

5.1  ARMS  -error-method  Statistics 

Table  XII  presents  the  overall  rms  errors  (both  initial  guess  and  final  pass)  for 
five  LRMW  parameters  for  twelve  observation  times  during  January  1963. 

Four  of  the  five  LRMW  parameter  verification  statistics  are  not  directly  com¬ 
parable  to  the  LMW  parameter  verification  statistics  because  three  of  them  are 
defined  differently  (as  discussed  in  Section  I)  and  thickness  of  the  LRMW  is  a  new 
parameter.  The  only  direct  comparison  with  the  LMW  verification  statistics  that  can 
be  made  is  for  the  wind  speed  maximum. 

Because  the  analysis  procedures  for  both  the  LRMW  and  the  LMW  are  identical 
for  the  wind  speed  maximum  after  the  initial-guess  field  has  been  generated,  the 
comparison  was  made  with  the  statistics  obtained  from  the  initial-guess  fields  (the 
January  1963  data  was  used  as  input  to  both).  Table  XIII  lists  the  analysis  and  with¬ 
held  station  rms  errors  for  the  initial  guess  of  Wg^  for  verification  Category  A 
(WgM  >  100  knots)  for  each  of  the  12  observation  times.  It  is  seen  that  the  LRMW 
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TABLE  XII 

OVERALL  RMS  ERRORS  FOR  LRMW  PARAMETERS 
FOR  12  OBSERVATION  TIMES 


Parameter 

Verification 

category 

Error  (initial  guess/final  pass) 

Total 

rms  error 

Analysis 
rms  error* 

Withheld 
rms  errort 

W«M 

A 

16.4 

16.5 

16.4 

SM 

13.4 

4.8 

18.4 

B 

12.4 

11.5 

13.2 

10.2 

3.0 

14.1 

C 

10.0 

9.7 

10.4 

7.2 

2.0 

10.0 

h(LR) 

A 

6021 

5440 

6552 

4523 

766 

6351 

B 

6126 

5499 

6695 

4146 

1225 

5734 

Z(LR) 

A 

4476 

4345 

4604 

3139 

774 

4372 

B 

4476 

3789 

5070 

3615 

790 

5051 

sb 

A 

2.88 

2.90 

2.85 

2.07 

0.66 

2.85 

B 

2.23 

2.19 

2.26 

1.61 

0.56 

2.21 

C 

1.52 

1.37 

1.66 

1.19 

0.31 

1.66 

s 

A 

2.81 

2.52 

3.08 

a 

2.14 

0.57 

2.97 

B 

1.91 

1.94 

1.87 

1.48 

0.47 

2.04 

C 

1.31 

1.26 

1.36 

1.03 

0.28 

1.43 

♦Verified back  to  stations  used  in  the  analysis. 

tVerified  to  those  stations  that  were  withheld  from  the  analysis. 
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equations  result  in  lower  rms  errors  for  11  of  12  times  for  the  analysis  stations  (in 
one  case,  the  error  is  the  same)  and  10  of  12  times  for  the  withheld  stations.  To  test 
the  significance  of  the  better  result  given  by  the  new  equations,  the  data  were  sub¬ 
jected  to  the  Student’s  t-test  for  paired  comparisons  [8].  It  was  found  to  be  significant 
at  less  than  the  1%  level  (i.e.,  highly  significant).  It  is  concluded  that  the  newly 
defined  jet-stream  categories  and  newly  derived  equations  produce  higher  quality 
initial  guesses  when  >  100  knots  than  the  LMW  equations.  For  the  lower  wind 

speed  verification  categories  (i.e.,  =  60—100  knots  and  <  60  knots),  the  differ¬ 

ence  in  the  LRMW  and  LMW  errors  for  the  initial-guess  fields  was  insignificant. 

For  the  other  four  LRMW  fields,  the  error  statistics  are  considered  to  be  low 
enough  to  indicate  that  the  analyses  are  of  high  quality  although  the  errors  for  h(LR) 
are  higher  than  one  might  hope  for.  The  computation  of  an  estimated  thickness  that 
may  produce  relatively  large  errors  over  about  20%  of  the  analysis  area  is  probably 
the  prime  reason  for  the  relatively  high  errors.  The  discussion  of  the  analyses  of  the 
LRMW  parameters  presented  in  the  next  section  shows  an  example  of  a  h(LR)  map  and 
suggests  a  method  for  detecting  areas  where  large  errors  are  likely,  and  a  solution 
to  the  problem. 

The  verification  statistics  for  all  the  LRMW  parameters  show  that  the  reduction 
of  the  total  rms  error  from  the  initial  guess  to  the  final  pass  is  due  almost  entirely 
to  the  reduction  of  the  rms  error  at  analysis  stations.  At  times,  the  withheld  station 
rms  error  is  somewhat  higher  for  the  final  pass  than  it  was  for  the  initial  guess. 
Similar  results  were  obtained  for  the  LMW  technique.  Logical  reasons  were  presented 
in  the  previous  report  [6]  and  are  repeated  here: 

“(a)  The  analysis  station  rms  error  essentially  indicates  how  well 
the  analysis  technique  fits  the  observations  (and  these  stations  may  be 
considered  as  minimum  error  points  over  the  analysis  area). 

“(b)  The  withheld  stations,  which  approximate  maximum  error  points 
do  not  contribute  to  the  analysis,  by  definition.  A  station  withheld  over  a 
sparse  data  area  may  be  the  only  station  in  the  area  (e.g.,  a  stationary 
weather  ship  in  the  ocean)  and  the  analysis  over  that  area  will  remain 
unchanged  from  the  initial  guess  to  the  final  pass. 

“(c)  Withheld  stations  over  a  medium  or  high  data-density  area  are  not 
allowed  to  influence  the  analysis.  The  surrounding  analysis  stations,  however, 
may  change  the  analysis  to  the  extent  that  the  point  at  which  there  is  a  with- 
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held  station  observation  may  actually  be  in  less  agreement  with  the  final 
analysis  than  it  was  with  the  initial -guess  analysis.” 

One  of  the  deficiencies  of  the  previous  LMW  technique  was  the  difficulty  of 
specifying  the  height  of  the  level  of  maximum  wind,  Z(L), because  sometimes  high 
wind  speeds  were  evident  over  a  thick  layer.  This  was  especially  true  when  the  maxi¬ 
ma  were  between  50  and  100  knots  with  the  height  of  the  LMW  varying  widely  over 
adjacent  grid  points  or  stations.  Thus,  the  rms  errors  for  the  initial  guess  and  final 
pass  of  Z(L)  for  verification  categories  A  and  B  were:  A  (  >  100  knots)— 4818  ft  and 
3533  ft,  respectively,  B  (50—100  knots)— 7378  ft  and  5738  ft,  respectively. 

The  mean  height  of  the  layer  of  maximum  wind,  Z(LR),  is  easier  to  specify 
and  more  meaningful,  particularly  for  verification  Category  B,  as  is  apparent  from 
comparison  of  the  Z(LR)  errors  with  the  Z(L)  errors. 

For  the  shears  below  and  above  the  LRMW,  the  rms  errors  are  about  the  same 
as  those  that  were  associated  with  the  LMW.  It  is  indicative  of  a  very  good  specifica¬ 
tion  of  the  LRMW  shear  field  because  this  result  was  achieved  in  spite  of  the  fact  that 
they  are  vector  shears  for  the  layer  5,000  ft  below  and  above  the  LRMW’  instead  of 
scalar  shears  for  the  layer  10,000  ft  below  and  above  the  LMW.  The  difference  in  the 

shear  criteria  (LRMW  and  LMW)  result  in  higher  average  values  for  shears  (by 

-3  -1  -3  -1 

2  knots  10  ft  )  and  generally  higher  standard  deviations  (by  1—2  knots  10  '  ft  ) 

for  the  LRMW  shears  over  the  LMW  shears  (see  Table  V  of  this  report  and  Table  IV 

of  [6]). 

5.2  Objective  (computer)  versus  Subjective  (man)  LRMW  Analysis  Comparison 

Figures  10—14  illustrate  objective  and  subjective  analyses  of  LRMW  parameters. 
In  general,  these  charts  and  others  reveal  that  there  are  only  minor  differences  in  the 
analyses  over  the  dense  data  areas  (e.g.,  the  U.S.),  but  differences  are  more  notice¬ 
able  over  the  sparse  data  areas  (e.g.,  Canada).  These  differences  are  attributed 
primarily  to  the  human  analyst  drawing  to  the  available  data  which  results  in  a 
tendency  for  maxima  and  minima  to  be  positioned  close  to  the  observing  stations 
when,  in  reality,  they  may  be  some  distance  away. 

The  objective  thickness  of  the  LRMW  clearly  indicates  some  unrealistically 
strong  gradients  over  small  areas  of  the  western  Pacific  [Fig.  11  (a)].  These  are 
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Fig.  10(a).  Objective  maximum  wind  speed,  1200Z,  3  January  1963. 
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Fig.  10(b).  Subjective  maximum  wind  speed,  1200Z,  3  January  1963. 
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Fig.  11(a).  Objective  thickness  of  LRMW,  1200Z,  3  January  1963. 
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Fig.  11(b).  Subjective  thickness  of  LRMW,  1200Z,  3  January  1963. 
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Jet  cores 


Fig.  12(a).  Objective  mean  height  of  LRMW,  1200Z,  3  January  1963. 


C 


56 


— 


— 


* 


Fig.  12(b).  Subjective  mean  height  of  LRMW,  1200Z,  3  January  1963.  * 
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Fig.  13(a).  Objective  vector  shear  below  the  LRMW,  1200Z,  3  January  1963. 
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Fig.  13(b).  Subjective  vector  shear  below  the  LRMW,  1200Z,  3  January  1963. 


59 


Fig.  14(a).  Objective  vector  shear  above  the  LRMW,  1200Z,  3  January  1963. 
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Fig.  14(b).  Subjective  vector  shear  above  the  LRMW,  1200Z,  3  January  1963. 
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TABLE  XIII 

LRMW  VERSUS  LMW  COMPARISON:  CATEGORY  A  (WgM  >  100  knots) 
INITIAL-GUESS  VERIFICATION  STATISTICS 


Date 

Analysis 

LRMW 

Station  rms  error 

LMW 

Withheld  station  rms  error 

LRMW 

LMW 

Jan 

1 

14.3 

22.3 

20.2 

23.4 

Jan 

2 

16.3 

22.1 

11.9 

17.3 

Jan 

3 

27.1 

27.1 

9.2 

19.7 

Jan 

4 

19.3 

22.4 

21.2 

20.2 

Jan 

5 

16.5 

24.4 

20.3 

26.3 

Jan 

6 

18.2 

21.2 

21.1 

26.3 

Jan 

7 

18.7 

28.0 

13.9 

15.6 

Jan 

8 

15.0 

30.0 

13.2 

23.6 

Jan 

9 

5.8 

7.6 

18.4 

19.8 

Jan 

12 

13.7 

27.0 

18.8 

20.7 

Jan 

13 

15.2 

21.4 

14.0 

11.0 

Jan 

14 

15.2 

19.7 

13.5 

19.6 

entirely  due  to  a  value  at  one  grid  point  that  is  in  considerable  variance  with  the 
surrounding  grid  points.  The  average  of  the  four  points  surrounding  the  central  point, 
for  each  of  the  four  areas  of  strong  gradients,  is  >  11,000  ft  different  than  the  value 
for  the  central  point.  It  is  very  likely  that  these  are  instances  where  the  wrong  category 
equation  is  used,  as  discussed  in  Sections  II  and  HI.  The  analyses  over  these  areas  can 
be  subjectively  (or  objectively)  adjusted  to  be  more  consistent  with  the  surrounding  grid 
points  and  a  high  degree  of  confidence  may  be  placed  in  the  adjusted  analysis  if  the 
strong  gradient  is  due  to  a  single  grid  point  value  (as  it  is  for  these  cases). 

5.3  Consistency  between  LRMW  Parameters 

A  study  of  the  LRMW  objective  analyses  showed  that  they  were  consistent  with 
one  another.  For  example:  (a)  the  maximum  vertical  shears  correspond  closely 
to  the  jet  cores  [Figs.  13  (a)  and  14  (a)];  (b)  the  mean  height  of  the  LRMW  is  at  low 
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elevation  for  the  mid-  and  high-latitude  polar  jets  (as  well  as  for  polar  jets  at  the 
lower  end  of  deep  troughs  that  extend  into  the  subtropics)  [Fig.  12(a)],and  (c)  the  mean 
height  of  the  LRMW  is  at  higher  elevations  for  subtropical  jets  in  the  lower  latitudes 
(as  well  as  for  jets  located  to  the  north  of  mid-  or  high-latitude  anticyclones)  [Fig. 
12(a)]. 

The  problem  of  low  shear  values  sometimes  appearing  near  jet  cores  in  the 
analysis  of  the  LMW  parameters  [6]  was  largely  eliminated  by  the  specification  of 
shears  below  and  above  the  LRMW  instead  of  below  and  above  the  LMW .  (The  LMW 
was  difficult  to  specify  in  cases  of  thick  vertical  layers  of  high  wind  speed,  and  either 
the  shear  below,  above,  or  both  were  sometimes  computed  to  be  relatively  low  values 
depending  on  where  in  the  layer  the  level  of  maximum  wind  was  specified.) 

The  thickness  of  the  LRMW  is  usually  <  10,000  ft  near  the  jet  cores  for  the  series 
of  maps,  and  thickness  maxima  are  often  located  between  two  jet  cores;  this  is  in 
agreement  with  what  is  known  to  occur.  For  example,  over  the  north  central  U.S.  on 
January  3,  1963  at  1200Z  there  is  a  polar  jet  core  [Z(LR)  ~  32,000  ft]  and  the  sub¬ 
tropical  jet  is  located  over  the  south-central  U.S.  and  extends  into  the  Central  Gulf  of 
Mexico  [Z(LR)  ~  40,000  ft].  The  thickness  of  the  LRMW  shows  minima  (~  8,000  ft) 
along  both  cores  and  a  maximum  of  16,000  ft  between  the  two  cores,  probably  indicating 
the  presence  of  both  jets  in  this  area,  with  the  Z(LR)  at  about  36,000  ft. 

On  the  wind  speed  maximum  analysis  charts,  the  generated  observations  along 
the  jet  cores  tend  to  elongate  the  maximum  wind  isotachs  from  those  isotachs  on  the 
initial-guess  chart.  The  purpose  of  the  generated  observations  is  to  give  a  better 
representation  of  the  wind  field  in  the  vicinity  of  the  jet.  In  some  areas  on  some 
charts  in  the  series,  withheld-station ..  observations  near  jet  cores  permitted  a  check 
on  whether  this  was  being  carried  out.  It  was  determined  that  for  most  cases  the 
generated  observations  did  give  a  better  representation  of  the  wind  field. 
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SECTION  IV 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  definition  and  objective  analysis  of  a  layer  of  maximum  wind  (LRMW)  and 
its  associated  parameters  of  wind  speed  maximum,  thickness,  mean  height,  and  vector 
shears  below  and  above  the  LRMW  leads  to  a  realistic  three-dimensional  picture  of 
the  jet  stream. 

The  categorization  procedures  for  wind  profiles  that  result  in  nine  jet-stream 
categories  for  which  initial-guess  regression  equations  were  derived  provide  high 
quality  initial-guess  fields  for  the  LRMW  parameters. 

A  significant  improvement  is  achieved  in  the  initial -guess  specification  of  the 
wind  speed  maximum,  in  areas  where  the  speed  exceeds  100  knots,  over  that  given  by 
the  level  of  maximum  wind  (LMW)  equations  [6J.  This  is  due  to  the  refinement  in 
the  categorization  procedure  that  includes  an  estimated  thickness  of  the  LRMW  as  part 
of  the  definition  criteria. 

Tests  of  the  LRMW  equations  on  independent  data  indicate  that  they  are  stable 
and  capable  of  specifying  very  well  the  general  characteristics  of  the  LRMW  profile. 

The  ability  of  the  analysis  technique  to  locate  jet  cores  and  generate  “observa¬ 
tions”  along  core  segments  results  in  better  definition  of  the  wind  field  in  the  vicinity 
of  the  jet  stream  than  that  given  by  the  initial-guess  “analysis”.  Over  the  data  areas, 
where  a  comparison  between  objective  and  subjective  analyses  can  be  made,  the 
objective  LRMW  analyses  compare  favorably  with  the  subjective  analyses. 

Because  the  LRMW  initial-guess  regression  equations  require  constant-pressure- 
surface  information  to  specify  the  LRMW  parameters,  they  are  purely  diagnostic 
equations,  but  may  be  used  with  prognostic  fields  to  produce  jet-stream  predictions. 

It  is  recommended  that  the  LRMW  diagnostic  equations  be  applied  to  a  series  of 
prognostic  constant-pressure-surface  fields  obtained  from  an  operational  numerical 
model.  The  resulting  LRMW  predictions  should  be  verified  and  evaluated.  If  they 
are  satisfactory,  the  equations  can  provide  a  relatively  simple  prediction  technique  for 
the  LRMW  parameters. 
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APPENDIX 


LRMW  INITIAL-GUESS  EQUATIONS 


1 


CATEGORY  1:  NO  JET 


W  =  12.865  +  0.39483 [W  (2)]  +  0.38694[W  (3)]  +  0.28719[W  (1)] 

SIVA  S  S  S 

S  =  1.4408  -  0. 0543 6 7 [W  (5-M)J  -  0.038540[W  (2-M)]  -  0.020353[W  (3-M)] 

0  s  s  s 

S’  =  0.17568  +  0.042254 [W  (1-M)]  +  0.054752[W  (2-M)]  +  0.026165[T(1)] 


CATEGORY  2:  MODERATE  SUBTROPICAL  JET 

W  =  9.8839  +  0.87219[Wg(2)]  +  0.15047[Ws(l)] 

h(LR)  =  -  2969.4  +  0.6818511^  +  3.8081[h  (3-2)] 

Z(LR)  =  4156.5  +  7. 6357 [W  (1)]  -  6.5161[W  (3)]  +  12.964[T(5)]  +  0.22620[h(LR)] 

s  s 

-  11.49  [Wg(2)]  +  9.2997[WgM] 

S’  =  12.569  -  0.11389[W  (3-M)]  -  0.0028825[Z(LR)]  -  0.036523[W  (5)] 

OS  s 

+  0.0010659  [h(LR)] 

Sa  =  8.2556  -  0.09569[WgM]  +  0.092069[Wg(l)]  -  0.0021269[Z(LR)] 


CATEGORY  3:  STRONG  SUBTROPICAL  JET 

WgM  =  17.385  +  0.90854[Wg(2)] 

h(LR)  =  1302.2  +  0.7  0960h^  +  18.102[T(2)] 

Z(LR)  =  3899.9  +  45.540[AT/An(2-1)]  -  4.7480[W  (3-M)]  +  0.19228[h(LR)] 

s 

+  7.9590[T(5)] 

Sb  =  -  1.0689  +  77.530S3  +  0.0019941[h(LR)]  -  0.037482[Wg(5-M)] 

-  0.02065[W  (3)] 
s 

S&  =  15.124  -  0.093949[WgM]  +  0.083590[Wg(l)]  -  0.00368[Z(LR)] 


CATEGORY  4:  SINGLE,  MODERATE  SUBTROPICAL-POLAR  JET 

W 


sM 


6.8307  +  0.49099[W  (2)]  +  0.55503[W  (3)]  +  0.15625[W  (1)] 
s  s  s 


h(LR)  =  706.54  +  16.776[W  (2-M)]  +  8.5505[W  (3)]  -  4.1051[W  (2)] 

s  s  s 

Z(LR)  =  4034.5  +  26.611[T(5)]  +  9.3408[W  (1)]  -  22.734[W  (3)]  +  20.991[W  (2)] 

s  s  s 

-  5.1075[WS(5)] 

S,  =  6.6065  -  0. 02205 [AW(5-3)]  -  0.0017313[Z(LR)]  +  0.094231[Wt 

-  0.10067[WS(5)] 

?a  =  -  0.48251  +  50.086S!  +  0.033009[WS(1-M)]  +  0.03936[Wg(3-M)] 


sMJ 
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CATEGORY  5:  SINGLE,  STRONG  SUBTROPICAL-POLAR  JET 


h(LR) 

Z(LR) 


S 

a 


12.373  +  0.67147[W  (3)]  +  0.42204[W  (2)] 
s  s 

529.61  +  0.40014h  +  4.2752[W  (2-M)] 

G  S 

-  1389.7  +  4.5992[h(5-3)]  +  l4.266[AW(3-2)]  -  0.35100[h(LR)] 
+  53 .791  [AT/  An(2-1)] 

-  0.36367  +  108. 92S  -  0.032949[W  (2-M)] 

o  s 

6.2663  +  77.157S  +  0.13921  [T(2)J  +  0.044702[W  (3-M)] 

1  s 


CATEGORY  6;  MULTIPLE,  MODERATE  SUBTROPICAL- POLAR  JE T 

WM  =  6.3829  +  0.50841[Wg(2)]  +  0.43791[Wg(3)]  +  0.11974[Wg(l)] 

h(LR)  =  1363.0  +  0.24537hj  +  3.1464[AW(7-5)]  +  2.5430[WT  (1-M)] 

G  S 

Z(LR)  =  4014.6  -  11.804 [W  (5)]  +  10.889[W  (1)]  +  16.757[T(5)]  -  23.391[W  (3)] 

s  s  s 

+  22.043 [W  (2)] 
s 

=  6.8636  +  126.47S5  +  0.0011678[h(LR)]  -  0.0012614[Z(LR)] 

-  0.02789[W  (7)]  +  0.053647[T(3)] 
s 

S’  =  -  0.85393  +  129. 06S, 
a  1 


CATEGORY  7:  MULTIPLE,  STRONG  SUBTROPICAL-POLAR  JET 


W  =  11.248  +  0.503 25 [W  (2)]  +  0.48457[W  (3)] 

SM  S  S 

h(LR)  =  -  1079.1  +  0.42786h  +  2.9815[W  (5)]  +  2.8140[AW(2-1)]  +  2.3720lh(3-2)] 

G  S 

Z(LR)  =  -  402.90  +  18.137[AW(3-2)]  +  4.1145[h(7-5)]  +  2.7292[W  (1)] 

s 

-  4 .3 7 02 [W  (5)]  +  0.29092[h(LR)]  -  13.770[W  (2-M)] 
s  s 

Su  =  9.7026  +  0. 084459 [A W(5-3)]  -  0.0023900[Z(LR)] 


0.29497  +  221. 09S, 


0.055363 [W  (1-M)]  -  0.015380[W  (5)] 
s  s 
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CATEGORY  8:  MODERATE  POLAR  JET 


W  w  =  3.4559  +  1.0407IW  (3)1 

sM  s'  ' 

h(LR)  =  176.37  +  0.57071h  +  4.9137[W  (5)] 

6  S 

Z(LR)  =  2988.4  +  25.356[T(5)]  -  8.8904[W  (5-M)]  +  4.0343[W  (1-M)] 

s  s 

-  9. 03 66 [T (2)] 

Sb  =  -  0.59366  +  111.33S5  +  0.0017857[h(LR)] 

=  0.71950  +  0.11020[Wg(2)]  -  0.094667WgM  -  0.0013345[Z(LR)] 

+  0.020708 [W  (1)] 
s 


CATEGORY  9:  STRONG  POLAR  JET 

=  9.7510  +  0.97260[W  (3)] 

s 

=  -  1635.3  +  0.64580h  +  2.9689[AW(3-2)]  +  1.8193[h(5-3)] 

=  3566.6  +  21.726[T(5)]  -  4.6704[W  (5)]  +  3.7669[W  (2)] 

s  s 

=  -  3.6220  +  152. 50S  +  0.0027131[h(LR)J 

o 

=  1.1766  +  0.12513 [W  (2-M)]  +  188.36S  -  0.064509[W  (1-M)] 

S  -L  S 


w  M 

sM 

h(LR) 

Z(LR) 
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